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Abstract 

Fractional classical mechanics has been introduced and developed 
as a classical counterpart of the fractional quantum mechanics. La- 
grange, Hamilton and Hamilton-Jacobi frameworks have been imple- 
mented for the fractional classical mechanics. The Lagrangian of frac- 
tional classical mechanics has been introduced, and equation of motion 
has been obtained. Fractional oscillator model has been launched and 
solved in ID case. A new equation for the period of oscillations of 
fractional classical oscillator has been found. The interplay between 
the energy dependency of the period of classical oscillations and the 
non-equidistant distribution of the energy levels for fractional quan- 
tum oscillator has been discussed. 

We discuss as well, the relationships between new equations of frac- 
tional classical mechanics and the well-known fundamental equations 
of classical mechanics. 

PACS numbers: 45.50.Dd; 45.05.+x; 45.20.Jj. 

Keywords: Fractional classical mechanics, Lagrange, Hamilton and 
Hamilton-Jacobi approaches, Fractional oscillator model, fractional 
Kepler's third law. 



1 Introduction 

We introduce and elaborate fractional classical mechanics as a classical coun- 
terpart of the fractional quantum mechanics To begin with, let us 
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consider the equation for canonical classical mechanics action Sa{p,q) (see, 
Eq.(23) in 0), 

tb 

Saip, q) = j dr{p{T)q{r) - if„(p(r), g(r))), (1) 

ta 

where the classical Hamiltonian Ha{p{T), q{T)) arrives from the Hamilto- 
nian Haip,q) 

H^ip,q) = D^\p\'' + Viq), 1< a < 2, (2) 

with V{q) as a potential energy and p — )■ p(r), q — )■ qi^r), and {p(t), g(r)} 
is the particle trajectory in phase space. The scaling factor Da has physical 
dimension [Da] = erg^~" ■ cm" ■ sec~", see, [I]- [2]. Since the Hamiltonian 
Haip, q) does not explicitly depend on the time, it represents a conserved 
quantity which is in fact the total energy of fractional classical mechanics 
system. 

Here we consider ID fractional classical mechanics. The 3D generalization 
is straightforward and it is based on the Hamiltonian [4] 

Ha{v>,ci) = Da\vr + V{ci), l<a<2, (3) 

where p and q are 3D vectors. 
The paper is organized as follows. 

Fundamentals of fractional classical mechanics are covered in Sec. 2, where 
Lagrange, Hamilton and Hamilton- Jacobi frameworks have been implemented. 

In Sec. 3 the scaling analysis of fractional classical motion equations has 
been developed based on the mechanical similarity. We discover and dis- 
cuss fractional Kepler's third law which is a generalization of the well-known 
Kepler's third law. 

The motion equations for the fractional classical mechanics have been 
integrated in Sec. 4. Fractional classical oscillator model has been introduced. 
We found an equation for the period of oscillations of fractional classical 
oscillator. The map between the energy dependence of the period of classical 
oscillations and the non-equidistant distribution of quantum energy levels 
has been established. 

In conclusion we outline the new developments. 
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2 Fundamentals of fractional classical mechan- 
ics 



2.1 The Lagrange approach 

The Lagrangian of fractional classical mechanics La{q, q) is defined as usual 

La{q,q) = pq - Ha{p,q), (4) 

where the momentum p is 

p=^^^, (5) 
dq 

and Ha{p, q) is given by Eq.(l2]). 

Hence, we obtain the Lagrangian of the fractional classical mechanics 

^.(g,9)= f^V"^l9l^-ng), 1<«<2. (6) 



For a free particle, V{q) = 0, the Lagrangian of the fractional classical 
mechanics is 

Further, the Euler-Lagrange equation of motion has the standard form 

d dLa{q,q) _ dLa{q,q) ^ ^ 
dt Qq dq 

By substituting Eq.Q into Eq.® we find the motion equation in La- 
grangian form 

1 1 -i-i^- dV(q) ^ 

+ o = 0- (9) 



^aDaJ a — 1 dq 

This equation has to be accompanied by the initial conditions. We impose 
the following initial conditions. At t = 0, the initial displacement is denoted 
by go and the corresponding velocity is denoted by qQ, that is. 
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q{t = 0)=qo, and q{t = 0) = q^. (10) 

As one can see, Eq.Q has a non- linear kinematic term. 
The new equation (Q is fractional generalization of the well known equa- 
tion of motion of classical mechanics in the Lagrange form. 
In the special case when a = 2, Eq.([9]) goes into 

^, + ^ = 0. (11) 
oq 

where m is a particle mass (-D2 = l/2m) and the initial conditions are 
given by Eq.( ITO|) . 

2.2 The Hamilton approach 

To obtain the Hamilton equations of motion for the fractional classical me- 
chanics we apply variational principle, 

6S^ip,q)=0, (12) 
where the action Sa{p, q) is given by Eq.([I]). 

Considering the momentum p and coordinate q as independent variables 
we can write 



SSM ,) = / dr(S„ + pS, - ^^^p - = 0. (13) 



Upon integration by parts of the second term p6q, the variation 6Sa 
becomes 



SS^ip, q)=j dr5p{q - ^^^) + p5qt - j dr5q{p + ^^^) = 0. 

Since 5q{ta) = Sq{tb) = at the end points of the trajectory, the term 
p6q is 0. Between the end points 6p and 6q can take on any arbitrary value. 
Hence, the variation 6Sa can be if the following conditions are satisfied 
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dp ^ dq ' 

This is, of course, the canonical Hamilton equations of motion. Thus, for 
the time- independent Hamiltonian given by Eq.(l2]) we obtain the equations 

q = aD^lpr'sgnp, p=-YM^ 1< a < 2, (15) 

where sgnp is the sign function which for nonzero values of p can be 
defined by the formula 

p 

sgnp = — , 
\P\ 

where \p\ is the absolute value of p. 

Equations ( !T5l) are the Hamilton equations for fractional classical me- 
chanics system. 

2.2.1 Poisson bracket approach to fractional classical mechanics 

It is well-known that Hamiltonian classical mechanics could be reformulated 
in terms of the Poisson brackets. For the two arbitrary functions u{p, q) and 
f (p, q) of variables p and the Poisson bracket is defined as [5] 

dudv dudv 
Wp,g),.(p,g)} = -----. (16) 

The Hamilton's equations of motion have an equivalent expression in 
terms of the Poisson bracket. Indeed, suppose that {{p,q,t) is a function of 
momentum p, coordinate q and time t. Then we have 

rff df fdf. m.\ 

dt = m^W^d-qV- ^'^^ 

Substituting q and p given by Eq.(lT4l) yields 

di (9f rrr n 



where {Ha,f} is the Poisson bracket. 
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{H^,n = ^--^-. (19) 

dp dq dq dp' 

with Ha given by Eq.(j2]). 
2.3 The Hamilton- Jacobi approach 

Having the Lagrangian ([6]), we can present classical mechanics action given 
by Eq.([T]) as a function of coordinate q, 



h h f 1 

SM) = j dTL^{q,q) = j dr\ (^-1-^" ^|g|^ - \/(g) j> , 1< a < 2, 

(20) 

Let's now treat the action Sa{q,h) as a function of coordinate q and 
the upper limit of integration, tb- Further, we suppose that motion goes 
along the actual path, that is we consider action as a functional of actual 
trajectory g(r) in coordinate space. To evaluate the variance 6Sa{q,tb) we 
have to compare the values of Sa{q, tb) for trajectories having a common start 
point at q(ta), but passing through different end points at the U- Thus, for 
the variation of the action by the variation of the end point Sq{tb) of the 
trajectory we have 



SSa{q,tb) = S I dTLa{q,q) = ^Sq\l + j dr5q f ^ - j (21) 

Because of Eq.® and Sq(ta) = we obtain 

5Saiq,tb) =p5q{tb), 

or 

= (22) 

For simplicity, letting 6q{tb) = 6q and tb = t, we can write 
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dSa{q,t) dSa{q,t) dSa{q,t) ■ dSa{q,t) . 

^^ = ^^ + ^^^=^^+^^- ^'^^ 

According to the definition of tlie action its total time derivative along 
the trajectory is dSa/dt = La, and Eq.(l23l) can be presented as 

^^^ = La-pq=-Ha{p,q), (24) 

where Eq.dl]) has been taken into account and p is given by Eq.(l22l). 
Therefore, we come to the equation 

^^ + ff<.(^^.,)=0. (25) 

at oq 

With Ha{p, q) given by Eq.(|2]) it takes the form 

5^ + Z,„|5«:M|<. + v(,)=0. l<a<2. (26) 

where 5*0, is called the Hamilton's principal function. 
For time-independent Hamiltonian the variables q and t in this equation 
can be separated. That is, we can search for solution of Eq. fl26|) in the form 



Sa{q,t,E) = S^^\q,E)-Et, (27) 

where Sa"* {q, E) is time-independent Hamilton's principal function and 
E, is the constant of integration, which has been identified with the total 
energy. Substituting Eq. (!27I) into Eq. (!26|) yields 

dS^^\q,E) l<a<2. (28) 

oq 

Equations fl26|) and ( 128|1 are the Hamilton- Jacobi equations of fractional 
classical mechanics. 

As an example, let us consider a free particle, V{q) = 0. The Hamilton- 
Jacobi equation Eq.( l26l) for a free particle is 

dSai,,t) ^ dSaiq,t) l<a<2. (29) 

at oq 

Therefore, from Eqs. fl27|) and fl28|) we obtain a solution to the Hamilton- 
Jacobi equation Eq 
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S^{q,t,E) = {^y/'^q-Et. (30) 

Follow the Hamilton- Jacobi fundamentals (see, for instance, [5]) we dif- 
ferentiate Eg. (pi]) over the energy E and put the derivative equal to a new 

constant 6, 
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dS^{q,t,E) 



— {—)^'\-t, (31) 



which yields 



and 



dE aD^^Da 



q = aD^{^y-^{t + S), (32) 



p=^^^^^|^ = (|-)^^ (33) 

where E > and 1 < a < 2. 

The equations (!32l) and ( l33l) define the trajectory of a free particle in 
fractional classical mechanics framework. 

In limit case when a = 2 and Da = l/2m, Eqs.©, (IIS!), ([26]) and ([28]) are 
transformed into the well-known Hamilton, Lagrange and Hamilton-Jacobi 
equations of classical mechanics for a particle with mass m moving in the 
potential field V{q). 



3 Mechanical similarity 

First of all, we intend to study the general properties of ID fractional classical 
motion without integrating motion equations ([9]), ([T5l) . and ([28]) . 

When the potential energy V{q) is a homogeneous function of coordinate 
q, it is possible to find some general similarity relationships. 

Let us carry out a transformation in which the coordinates are changed 
by a factor p and the time by a factor r: q ^ q' = pq, t ^ t' = rt. Then all 
velocities q = dq/dt are changed by a factor p/r, and the kinetic energy by 
a factor (p/r)"^*-""^-*. If V{q) is a homogeneous function of degree /3 then it 
satisfies, 
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V{pq) = p^V{q). (34) 

It is obvious that if p and r are such that (p/r)"^*'"""^^ = p^ , i.e. r = 
, then the transformation leaves the motion equation unaltered. This 
invariance is called mechanical similarity. 

A change of all the coordinates of the particles by the same scale factor 
corresponds to replacement of the classical mechanical trajectories of the par- 
ticles by other trajectories, geometrically similar but different in size. Thus, 
we conclude that, if the potential energy of the system is a homogeneous 
function of degree /3 in Cartesian coordinates, the fractional equations of 
motion permit a series of geometrically similar trajectories, and the times of 
the motion between corresponding points are in the ratio 

where j is the ratio of linear scales of the two paths. Not only the times 
but also any other mechanical quantities are in a ratio which is a power of 
y. For example, the velocities and energies follow the scaling laws 



and 



q fl'Y'^ E' _ //' 



t' E'^ 



(36) 



t \ E 



(37) 



The following are some examples of the foregoing. 

It follows from Eq.(|3711 that at ^ + = 1 the period of oscillations does 
not depend on energy of oscillator. The condition ^ + = 1 with 1 < a < 2, 
1 < /3 < 2 brings a = 2 and /3 = 2 only. It means, for instance, that 
considering the fractional classical oscillator model with Hamiltonian (HTj) . 
we can conclude that the standard classical harmonic oscillator only, 

^2(p,g)=pV2m + (?V, (38) 

has the period of oscillations which does not depend on energy of oscilla- 
tor. 
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In the uniform field of force, the potential energy is a linear function of 

the coordinates, i.e. (3 = 1. From Eq. fl37|) we have j = (j)^^"'- Therefore, 
the time of motion in a uniform field (/3 = 1) is as the a-root of the initial 
altitude. 

If the potential energy is inversely proportional to the distance apart, i.e. 
it is a homogeneous function of degree (3 = —1, then Eq.( l35|l becomes 



t \ I 



(39) 



For instance, regarding the problem of orbital motion in 3D fractional 
classical mechanics, we can state that the orbital period to the power a is 
proportional to the power 2q; — 1 of its orbit scale. This statement is in fact 
a generalization of the well-known Kepler's third law. We call Eq.( l39l) as 
fractional Kepler's third law. In special case a = 2, Eq. fl39|) turns into the 
well-known Kepler's third law [5]. 

In general, for negative (3, P = —7, in the space of the fractional Hamil- 
tonians Ha--y at 

^ + l-- = i (40) 
a L 

there exists the subset of fractional dynamic systems, orbital motion of 
which follows the well-known Kepler's third law, that is, the square of the 
orbital period is proportional to the cube of the orbit size. 

4 Integration of the motion equations 

4.1 Fractional classical ID oscillator 
4.1.1 The Lagrange approach 

Fractional quantum oscillator model has been introduced in pj. Fractional 
classical oscillator can be considered as a classical counterpart of the frac- 
tional quantum oscillator model. 

We introduce fractional classical ID oscillator as a mechanical system 
with the Hamiltonian 

ff(p,g) = Z},|pr + /|g|^ (41) 
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where g is a constant with physical dimension [g] = erg^/^ ■ cm~^^'^ and a 
and /3 are parameters, l<a<2, l</3<2. 

It follows from Eq.(j6]) that the Lagrangian of the fractional classical ID 
oscillator is 



1 

a — 1 



Liq,q)=[—-) l<a<2, 1< /3 < 2. (42) 

\aUaJ a 

Hence, the motion equation of fractional classical ID oscillator is 

1 

-g|g|i^ + /3/|g|'^-^sgng = 0. (43) 

This is a new non-linear classical mechanics equation of motion. When 
a = 2 and /3 = 2 it goes into the well-known linear equation of motion for 
classical mechanics oscillator which has the form, 

mq + 2g'^q = 0. (44) 
where m is a particle mass {D2 = l/2m). 

To integrate Eq. fH3l) we start from the law of energy conservation. For 
the Lagrangian Eq. (H2|) we have 

(^) "'''"^^ '''''' = '''' 

where E is the total energy. Thus, we have 

\q\=aDl/'^{E-g'\qf)^, (46) 

which is a first-order differential equation, and it can be integrated. 

Since the kinetic energy is positive, the total energy E always exceeds the 
potential energy, that is E > g'^\q\'^. The points where the potential energy 
equals the total energy g'^\q\'^ = E are turning points of classical trajectory. 
The ID motion bounded by two turning points is oscillatory, the particle 
moves repeatedly between those two points. 

Substituting q = {E/g'^Y^^y into Eq.l^ yields 



\y\ 



^^i/"^2//3^i-{i+i)(^ - \yfy-^. (47) 
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Hence, the period T(a, (3) of oscillations is 




By substituting z = we rewrite the last equation in the form 




(48) 



With the help of the 5-function definition [6] 
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5(i -) = [ dzz^-\l-z) 



1-1 





we finally find for the period of oscillations of fractional classical ID os- 
cillator 



tional classical ID oscillator. The dependency on energy is in agreement with 
the scaling law given by Eq. fl?r|) . 

It follows from Eq.( H9|) that period T(a, (3) doesn't depend on the energy 
of oscillator when (^ + -i) = l. Ifl<Q;<2 and 1 < (3 < 2 then condition 

+ ^) = 1 gives us that a = 2 and (3 = 2. Hence, we come to the 
standard classical mechanics harmonic oscillator with the Hamiltonian given 
by Eq. fl38|) and the energy independent oscillation period, T(2, 2) = n\/2m/g. 

Table 1 shows that the energy dependency of the period of fractional clas- 
sical oscillator is a classical counterpart of the fractional quantum mechanics 
statement on non-equidistant energy levels of fractional quantum oscillator. 
For classical mechanics, independence on energy of the period of classical 
oscillator is a classical counterpart of the quantum mechanics statement on 
equidistant energy levels of quantum oscillator. 




(49) 
;y of frac- 
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Period of oscillations T 


Energy levels Sn 


l<a<2 
l</3<2 


Fractional classical mechanics 


i<'ractional quantum mechanics 
[ 2S(^ 1 + 1) 1 ^"+2) 


a=2 
/3=2 


Classical mechanics 
n\/2m/g 


Quantum mechanics 



Table 1. Energy dependence/independence of the period of classical oscil- 
lator vs non- equidistant /equidistant energy levels of quantum oscillato-^. 



4.1.2 The Hamilton approach 

For the fractional ID oscillator the Hamilton equations of motion in accor- 
dance with Eqs.( !T4|) and (14T]) are 



dp 



aDM^-hgnp, (50) 



P = -^ = -/3/|gr^sgng, (51) 



dq 

where l<a<2, l</3<2. 

Further, Hamilton equation ( 150|1 leads to 



Q-l 



l^l=(^j I'l^- (52' 

Substituting Eq. fl52l) into Eq. fj^T]) yields exactly Eq.f HSj) . Therefore, down- 
stream integration can be done by the same way as it was done above in the 
framework of the Lagrange approach. 



4.1.3 The Hamilton- Jacobi approach 

For fractional classical oscillator with the Hamilton function defined by Eq. 0411) 
a complete integral of the Hamilton- Jacobi equation fl27|) is 

"'^The period of oscillations T and frequency of oscillations uj are related by T = ^ixjuj. 

Hence, the classical mechanics frequency is w = ^^^9^ ^'^'^ terms of the frequency, the 

quantum mechanic energy levels are En — fioj{n + 

The energy levels equation for the fractional quantum oscillator has been taken from 

Il- 



ls 



S{q,t,E) = I dqi^jE - g'\qf)y/^ - Et, (53) 

where the integration constant E can be identified with the total energy 
of the fractional classical oscillator. 

It is known that in standard classical mechanics [5], a = 2, the integral 
in Eq.( l53|) is considered as an indefinite integral. At this point, to move 
forward, we introduce the ansatz to treat the integral in Eq. (l53|) as 

S{q,t,E) = j dq{^{E - g^\q\^)Y'^ - Et. (54) 







Following the Hamilton- Jacobi fundamentals (see, for instance, ^) we 
differentiate Eq. (IS^ over the energy E and put the derivative equal to a new 
constant 5, 



dS{q,t,E) _ 1 
dE " 



t^g'(E-/|g'|^)(i-")/"-t = 5. (55) 







Substituting integration variable q' with new variable z, q' = {E / g'^y/^ z^^^ 
yields 

t + ^=—-U7r-r-. / dzz-^-\l~z)^'\ (56) 



apoTg'"'^ 

The integral in the above equation can be expressed as the incomplete 
Beta functiorl^ and we obtain 

(- + -)-! 

'''' 

This equation is the solution of the fractional ID oscillator as function 
t{q) in terms of the incomplete Beta function. 



^The incomplete Beta function is defined as [7] 

X 

B.ifi,!.)^ J dvyf^-^l ~ yr-\ (57) 
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The incomplete Beta function Bx{^, v) has the hypergeometric represen- 
tation 

B^{lj,,u) = — 1 - z/;/i + l;a;), (59) 
A'' 

where F{fi, 1 — u; fi + 1; x) is the hypergeometric function [8]. 
Now we can write Eq.( l58|) in the form 

Thus, in terms of the hypergeometric function we found solution of the 
fractional ID oscillator as function t{q). 

Let us show that the ansatz given by Eg. (15^ and based on it solution 
given by Eq. (!60!) allow us to reproduce the well-known solution to standard 
classical harmonic oscillator, a = 2 and /3 = 2. Indeed, at a = 2 and /3 = 2 
we have 

io{t + 6)=xF{^l^-;x% (61) 
where a new variable x has been introduced as, 

X = g^/l/Eq, (62) 

and uo = ^j2/mg is the frequency of classical ID oscillator. Since xF{}^, i; | 
arcsinx, it follows from Eq. (l6T]) that 

x{t) = sinu;(t + 6). 
By restoring the original dynamic variable q{t) from Eq.f p^ we obtain. 



2E 

(lit) = \l ^sincj(t + 5), (63) 



which is the solution of the Hamiltron-Jacobi equation for standard ID 
harmonic oscillator with the Hamiltonian given by Eq.( l38|) . 

Thus, we proved that the ansatz given by Eq. ( 15^ gives us the well-known 
solution to standard classical mechanics harmonic oscillator. 
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5 Conclusion 



Fractional classical mechanics has been introduced as a classical counterpart 
of fractional quantum mechanics. The Lagrange, Hamilton and Hamilton- 
Jacobi frameworks have been developed for fractional classical mechanics. 
Scaling analysis of fractional classical motion equations has been imple- 
mented based on the mechanical similarity. We discover and discuss frac- 
tional Kepler's third law which is a generalization of the well-known Kepler's 
third law. 

Fractional classical oscillator model has been introduced and motion equa- 
tions for the fractional classical oscillator have been integrated. We found an 
equation for the period of oscillations of fractional classical oscillator. The 
map between the energy dependence of the period of classical oscillations and 
the non-equidistant distribution of the energy levels for fractional quantum 
oscillator has been established. 

In the case when a = 2, all new developments are turned into the well- 
known results of the classical mechanics. 
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